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The classical experiments on turbulent friction in rough pipes were performed by J. Nikuradse in 
the 1930's. Seventy years later, they continue to defy theory. Here we model Nikuradse's experiments 
using the phenomenological theory of Kolmogorov, a theory that is widely thought to be applicable 
only to highly idealized flows. Our results include both the empirical scalings of Blasius and Strickler, 
and are otherwise in minute qualitative agreement with the experiments; they suggest that the 
phenomenological theory may be relevant to other flows of practical interest; and they unveil the 
existence of close ties between two milestones of experimental and theoretical turbulence. 



Turbulence is the unrest that spontaneously takes over 
a streamline flow adjacent to a wall or obstacle when 
the flow is made sufficiently fast. Although most of the 
flows that surround us in everyday life and in nature 
are turbulent flows over rough walls, these flows have 
remained amongst the least understood phenomena of 
classical physics Thus, one of the weightier ex- 

perimental studies of turbulent flows on rough walls, and 
the most useful in common applications, is yet to be ex- 
plained theoretically 70 years after its publication. In 
that study Nikuradse elucidated how the friction co- 
efficient between the wall of a pipe and the turbulent 
flow inside depends on the Reynolds number of the flow 
and the roughness of the wall. The friction coefficient, 
/, is a measure of the shear stress (or shear force per 
unit area) that the turbulent flow exerts on the wall of 
a pipe; it is customarily expressed in dimensionless form 
as / = r/yoF^, where p is the density of the liquid that 
flows in the pipe and V the mean velocity of the flow. 
The Reynolds number is defined as Re = VR/v, where 
R is the radius of the pipe and u the kinematic viscosity 
of the liquid. Last, the roughness is defined as the ratio 
r/R between the size r of the roughness elements (sand 
grains in the case of Nikuradse's experiments) that line 
the wall of the pipe and the radius of the pipe. 

Nikuradse presented his data in the form of six curves, 
the log-log plots of / versus Re for six values of the rough- 
ness These curves are shown in Fig.^ At the onset 
of turbulence 0, at a Re of about 3,000, all six curves 
rise united in a single bundle. At a Re of about 3,500, 
the bundle bends downward to form a marked hump and 
then it plunges in accord with Blasius's empirical scal- 
ing [3i / ~ Re"^^''', as one by one in order of decreas- 
ing roughness the curves start to careen away from the 
bundle. After leaving the bundle, which continues to 
plunge, each curve sets out to trace a belly @] as it steers 
farther from the bundle with increasing Re, then flexes 
towards a terminal, constant value of / that is in keep- 
ing with Strickler's empirical scaling j], / ~ {r/RY^^. 
For seventy years now, our understanding of these curves 
has been aided by little beyond a pictorial narrative of 




FIG. 1: Nikuradse's data. Up to a Re of about 3,000 the 
flow is streamline (free from turbulence) and / ~ 1/Re. Note 
that for very rough pipes (small R/r) the curves do not form 
a belly at intermediate values of Re. Inset: veriflcation of 
Strickler's empirical scaling for / at high Re, / ~ {r/R)^''^. 



roughness elements being progressively exposed to the 
turbulent flow as Re increases (sj. 

In our theoretical work, we adopt the phenomeno- 
logical imagery of "turbulent eddies" 0, 0, 0| and 
use the spectrum of turbulent energy at a length- 
scale a, E{a), to determine the velocity of the eddies 
of size s. Us, in the form = E{a)(j^'^da, where 
E{a) = Ae'^/^a^/^Cd{v/cr)ce{(7/R)- Here A is a di- 
mensionless constant, e is the turbulent power per unit 
mass, rj = i/3/4£-i/4 jg ^^iq viscous lengthscale, R is the 
largest lengthscale in the flow, Ae^/'^cr^/'^ is the Kol- 
mogorov spectrum (which is valid in the inertial range, 
rj <^ a <^ R), and Cd and Cg are dimensionless corrections 
for the dissipative range and the energetic range, respec- 
tively. For Cd we adopt an exponential form, Cdirj/a) — 
e:xjp{— Pt] / a) (which gives f» 1 except in the dissipa- 
tive range, where a ^ if), and for Ce the form proposed 
by von Karman, Ce{a / R) = (1 7(cr/i?)2)-iV6 (which 
gives Ce ~ 1 except in the energetic range, where a ^ R), 
where /3 and 7 are dimensionless constants To ob- 
tain expressions for Us and 77 in terms of Re, r/R, and 
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V, we invoke the usual scalings fis'l, e = k^u^/R (Tay- 
lor's scaling 0|, where is the characteristic velocity of 
the largest eddies and a dimensionless constant) and 
ur = (where dimensionless constant). Then, 

we can write rj = 6i?Re~^^^, where b = {KgK^)~^/'^, and 
(after changing the integration variable to a; = a/R) 
= AkJ\IJ^^'^ x'^/^CdibRe-"^/^ /x)ce{x)dx. For 
s ^ R we can set Ce — I, compute the integral, and 
let Re ^ oo to obtain ^ (3/2)A k^^^^ s/i?)2/3, 
or ul ^ uj^{s/R)'^/^, a well-known result of the phe- 
nomenological theory. Further, for consistency with 
Taylor's scaling we must have Ak^^^ = 2/3 (so 
that Ms — Ur for s ^ R) and therefore uj. — 
kI y2(2/3)/;/^x-i/3cd(foRe-3/Vx)ce(x)dx. 

We now seek to derive an expression for r, the shear 
stress on the wall of the pipe. We assume a viscous layer 
of constant thickness ary, where a is a dimensionless con- 
stant, and call W a wetted surface parallel to the peaks 
of the viscous layer (Fig. [^J. Then, r is effected by mo- 
mentum transfer across W. Above W, the velocity of the 
flow scales with and the fluid carries a high horizon- 
tal momentum per unit volume (~ pV). Below W, the 
velocity of the flow is negligible, and the fluid carries a 
negligible horizontal momentum per unit volume. Now 
consider an eddy that straddles the wetted surface W. 
This eddy transfers fluid of high horizontal momentum 
downwards across W, and fluid of negligible horizontal 
momentum upwards across W. The net rate of transfer 
of momentum across W is set by the velocity normal to 
W, which velocity is provided by the eddy. Therefore, 
if Vn denotes the velocity normal to W provided by the 
dominant eddy that straddles W, then the shear stress 
effected by momentum transfer across W scales in the 
form T ^ p Vvn ■ 

In order to identify the dominant eddy that straddles 
W, let us denote hy s = r + arj the size of the largest 
eddy that flts the coves between successive roughness el- 
ements. Eddies much larger than s can provide only a 
negligible velocity normal to W . (This observation is 
purely a matter of geometry.) On the other hand, eddies 
smaller than s can provide a sizable velocity normal to 
W . Nevertheless, if these eddies are much smaller than 
s, their velocities are overshadowed by the velocity of the 
eddy of size s. Thus, u„ scales with u^, which is the ve- 
locity of the eddy of size s, and the dominant eddy is 
the largest eddy that flts the coves between successive 
roughness elements. We conclude that r ~ pVus^ or 
T = KtP Vug (where Kt is a dimensionless constant of 
order 1), and therefore / — KtUs/V or 

/ rs/R \ 
f = K[j x-'/^CdibRe-^/^/x)ceix)dx\ , (1) 

where K = KrKu\/V^, s/R = r / R + abRe~^'^ , and b = 
(KeK^J^i/"*. Equation gives / as an explicit function 




FIG. 2: Schematic of the immediate vicinity of the wall with 
roughness elements of size r covered by a viscous layer of uni- 
form thickness arj. The distance between roughness elements 
is about equal to the hei ght of the roughness elements, as in 
Nikuradse's experiments |2C| . The horizontal line is the trace 
of a wetted surface W tangent to the peaks of the viscous 
layer. 

of the Reynolds number Re and the roughness r/R. 

To evaluate computationally the integral of 1^, we set 
(3 = 2.1, 7 = 6.783 (the values given in lU]), a = 5 
(brj being a common estimation of the thickness of the 
viscous layer), = 5/4 (a value that follows from Kol- 
mogorov's four-fifth law |ia]), k„ = 0.036 (0.036 ± 0.005 
being the value measured in pipe flow by Antonia and 
Pearson li*!), b = {k^k^)~^/'^ = 11.4, and treat Kt as a 
free parameter (albeit a parameter constrained by the- 
ory to be of order 1). With Kt = 0.5 (and therefore 
K = 0.015), (d gives the plots of Fig. |3 (Note that a 
different value of Kr would give the same plots except 
for a vertical translation.) These plots show that is 
in excellent qualitative agreement with Nikuradse's data, 
right from the onset of turbulence, including the hump 
and, for relatively low roughness, the bellies. These plots 
remain qualitatively the same even if the value of any of 
the parameters is changed widely. In particular, there 
is always a hump and there are always bellies: these are 
robust features closely connected with the overall form of 
the spectrum of turbulent energy. The connections will 
become apparent after the discussion that follows. 

To help interpreting our results, we compute / without 
including the correction for the energetic range — that is, 
setting 7 = 0. In this case, the integral of may be 
evaluated analytically, with the result 

/ = X(r/i? + a6 Re-3/4)i/3y;p(^^ (2) 

where F{y) — y'^^^^-2/3{y): ^-2/3 is the gamma func- 
tion of order -2/3, and y = (3t]/s = /36Re^^/''(r/i? + 
a6Re~"^/'*)^i. With the same values of Kt, k„, a, b, and 
(3 as before, 10) gives the solid-line plots in the inset of 
Fig. 13 The hump is no more. We conclude that the 
hump relates to the energetic range. Further, with the 
exception of the hump at relatively low Re, the plots of 
coincide with the plots of ((2Jl; thus, we can study 
to reach conclusions about Q at intermediate and high 
Re. For example, Q gives / ~ {r/R)'^/'^ for r '3> arj and 
/ - Re"^/^ for r < 077. It follows that both Q and ((TJ 
give a gradual transition between the empirical scalings 
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FIG. 3: Plot of Q. Inset: Plot of ^ (no correction for the 
energetic range: solid lines) and plot of JJJ with 7 = (no 
correction for the energetic range and the dissipative range: 
dashed lines). 



of Blasius and Strickler [l3|, in accord with Nikuradse's 
data. 

If we set /? = in addition to 7 = 0, (O simplifies to 
/ = KrKu {r/R + a^Re^^/*)^/^. With the same values 
of Kt, k„, a, and 6 as before, this expression gives the 
dashed-line plots in the inset of Fig. 13 Now the bellies 
are no more. We conclude that the bellies relate to the 
dissipative range. The dissipation depresses the values of 
/ at relatively low and intermediate Re, leading to the 
formation of the bellies of Nikuradse's data. 

We are ready to explain the unfolding of Nikuradse's 
data in terms of the varying habits of momentum trans- 
fer with increasing Re (Fig.0J. At relatively low Re, the 
inertial range is immature, and the momentum transfer is 
dominated by eddies in the energetic range, whose veloc- 
ity scales with V , and therefore with Re. Consequently, 
an increase in Re leads to a more vigorous momentum 
transfer — and to an increase in /. This effect explains 
the rising part of the hump. At higher Re, the momen- 
tum transfer is dominated by eddies of size s ~ ary 3> r. 
bmce rj ~ Re-^/^ with increasing Re the momentum 
transfer is effected by ever smaller (and slower) eddies, 
and / lessens as Re continues to increase. This effect ex- 
plains the plunging part of the hump — the part governed 
by Blasius's scaling. At intermediate Re, s — r^a-q with 
r « arj. Due to the decrease in rj, s continues to lessen as 
Re continues to increase, but at a lower rate than before, 
when it was s arj ^ r. Thus, the curve associated 
with r deviates from Blasius's scaling and starts to trace 
a belly. As 77 continues to decrease, the dominant ed- 
dies become decidedly larger than the smaller eddies in 
the inertial range, which is well established now, and any 
lingering dissipation at lengthscales larger than s must 
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FIG. 4: Schematic of the relations among a generic Niku- 
radse curve, the spectrum of turbulent energy, the size of the 
roughness elements, the thickness of the viscous layer, and 
the size of the dominant eddies. 



cease. This effect explains the rising part of the belly. 
Last, at high Ke, s ^ r ^ arj. As Re increases fur- 
ther, 1] lessens and new, smaller eddies populate the flow 
and become jumbled with the preexisting eddies. Yet 
the momentum transfer continues to be dominated by 
eddies of size r, and / remains invariant. This effect ex- 
plains Nikuradse's data at high Re, where / is governed 
by Strickler's scaling. 

We have predicated equation on the assump- 

tion that the turbulent eddies are governed by the phe- 
nomenological theory of turbulence. The theory was 
originally derived for isotropic and homogeneous flows, 
but recent research ISj suggests that it applies to much 
more general flows as well. Our results indicate that even 
where the flow is anisotropic and inhomogeneous — as is 
the case in the vicinity of a wall — the theory gives an ap- 
proximate solution that embodies the essential structure 
of the complete solution (including the correct scalings of 
Blasius and Strickler) and is in detailed qualitative agree- 
ment with the observed phenomenology. Remarkably, 
the qualitative agreement holds starting at the very onset 
of turbulence, in accord with experimental evidence that 
"in pipes, turbulence sets in suddenly and fully, without 
intermediate states and without a clear stability bound- 
ary" . The deficiencies in quantitative agreement point 
to a need for corrections to account for the effect of the 
roughness elements on the dissipative range as well as for 
the effect of the overall geometry on the energetic range. 

In conclusion, to a good approximation the eddies in 
a pipe are governed by the spectrum of turbulent energy 
of the phenomenological theory. The size of the eddies 
that dominate the momentum transfer close to the wall 
is set by a combination of the size of the roughness ele- 
ments and the viscous lengthscale. As a result, the de- 
pendence of the turbulent friction on the roughness and 
the Reynolds number is a direct manifestation of the dis- 
tribution of turbulent energy given by the phenomeno- 
logical theory. This close relation between the turbu- 
lent friction and the phenomenological theory |l9l | may 
be summarized in the following observation: the simi- 
larity exponents of Blasius and Strickler are but recast 
forms of the exponent 5/3 of the Kolmogorov spectrum. 
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